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Black hole evaporation may lead to massive or massless remnants, or naked singularities. 
This paper investigates this process in the context of two quite different two dimensional 
black hole models. The first is the original CGHS model, the second is another two dimen- 
sional dilaton-gravity model, but with properties much closer to physics in the real, four 
dimensional, world. Numerical simulations are performed of the formation and subsequent 
evaporation of black holes and the results are found to agree qualitatively with the exactly 
solved modified CGHS models, namely that the semiclassical approximation breaks down 
just before a naked singularity appears. 
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1. Introduction 

Ever since Hawking's discovery of tlie thermal cliaracter of radiation from black 
holes, much effort has been expended trying to understand the implications of this for 
a theory of quantum gravity. Conservative suggestions have been that uncharged black 
holes evolve into massive states somehow made stable by quantum gravity corrections, 
or into massless remnants such as the cornucopions, in which information is stored in a 
semi- infinite horn geometry. Both these possiblities allow unitarity to be preserved 0. 
A radical suggestion was made by Hawking to describe the physics of the endpoint of 
black hole evaporation. It was suggested that pure quantum states must evolve into mixed 
states, leading to a breakdown in predictability. This is an inevitable consequence of 
the formation of a naked singularity when the event horizon ceases to exist as the black 
hole evaporates away to leave empty space behind. This of course violates the Cosmic 
Censorship hypothesis of Penrose and leaves us unable to predict the future even at the 
classical level in half of spacetime. 

Recently much progress has been made in developing toy models of black hole evap- 
oration in two dimensions. This was initiated by the work of CGHS 0]. The quantum 
backreaction is included at the one loop level by integrating out the matter fields, in- 
ducing a Polyakov-type term in the effective action. The large limit is then taken 
allowing the one loop effects of the gravitational degrees of freedom to be ignored. Static 
solutions of the resulting equations have been obtained in f^HH- More recently analytic 
results for dynamical shock wave solutions have been obtained in modified models based on 
a Liouville type field theory [1T0| - 12 1 . There it was found that a naked singularity appears 
as the black hole continues evaporating to negative mass. 

This paper begins by reviewing the results of CGHS and then numerical results are 
presented for shock wave solutions in the original CGHS model. They are found to be in 
qualitative agreement with the analytic results of the modified models, namely that the 
semiclassical approximation breaks down just before a naked singularity appears. 

A model is then introduced which is obtained by the dimensional reduction of the four 



dimensional Einstein action. See for a discussion of a general class of two dimensional 
dilaton-gravity theories. This model is then coupled to matter generally covariant only in 
two dimensions allowing the backreaction to be included in the same way as the CGHS 
model. The black hole solutions are of course essentially the same as their four dimensional 
analogs, so it is hoped that the study of this model will lead to greater insights into physics 
in the real world. In particular the Hawking temperature is T = I/SttM. 
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This model seems even more difficult than the CGHS model to treat analytically, so 
once again numerical methods are used to study both the static solutions and dynami- 
cal shock waves. The finite mass static solutions exhibit singular horizons, reminiscent 
of the solutions of the CGHS model [|^,|^ . Again shock wave solutions evolve until a 
naked singularity appears, at which point the approximation breaks down, and fiuctua- 
tions in the gravitational degrees of freedom can not be neglected. At the very least this 
means observers may see high energy quantum gravity effects at the endpoint of black hole 
evaporation, but the question of what really happens is unresolved. 



2. Review of CGHS model 

CGHS begin with the classical action for dilaton gravity in two spacetime dimen- 
sions 

1 r [ 1 ^ 

Sd = - (fa e-^"^ -2a+a_p + 4a+(/>a_(/>-AV (2.1) 

where (f), p and / are the dilaton, Liouville and matter fields respectively. 

Static solutions of this action representing two dimensional black holes are 

e-2<^ = e-2^ = ^-A2x+x-. (2.2) 
A 

Here M corresponds to the ADM mass of the solution. The case M = corresponds to a 
coordinate transformation of the linear dilaton vacuum solution 

p = and (f) = — Aa , (2-3) 

where a = |(cr+ — a~) and 

e^^^ = Ax+, e"^^" = -Ax" . (2.4) 

Shock waves corresponding to a black hole with mass M being formed by collapse take the 
form 

g-2</> _ g-2p _ -x^x+x- - ^(a;+ - x+)e{x-^ - x+) . (2.5) 

\Xq 

The key point of 0] is that in the one loop approximation quantum corrections to 
the classical action may be computed simply by integrating out the fiuctuations in the 
matter fields. Taking the large limit, where A^ is the number of matter fields means the 
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contributions of the gravitational degrees of freedom may be ignored. The resulting action 
is 

Seff = Sci + lj d'^a ^d+pd.p . (2.6) 
The equations of motion for the case fi=0 can be written in the convenient form 



(2.7) 



It should be noted ( |2.7| ) become degenerate when cj) = (per = —\ log(A^/12), where 
in general a curvature singularity will appear. However, one expects the large approxi- 
mation to break down before this singularity is reached. This may be seen by examining 
the matrix K of the dilaton-Liouville kinetic term K{(f)) c)_$ where $ denotes the 
2- vector (p, (f)). The determinant of this matrix is 

deti^ = e-2'^(e-2'^ - ^) . (2.8) 

This should be 0{N'^) if the one loop contributions of the dilaton and Liouville fields are to 
be ignored, so one sees the large approximation fails when (f)cr — 4'^0{1). The quantity 
^2^"^^ is therefore a useful measure of the coupling of quantum gravity. When it becomes 
of order one, quantum gravity effects become strong. 

The following formula for the effective mass of a shock wave solution will be useful in 
the next section 

Me// = ^(l-^e2^)^e-2^i?. (2.9) 
This agrees precisely with the Bondi mass of a shock solution along the infall line a;"*" 



X 
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even in regions of strong coupling. Also in the limit x~ — oo this agrees with the 
total initial energy of the matter shock wave. After formation, i.e. for > xj, M^ff 
will decrease as x~ increases corresponding to energy loss as the black hole evaporates. 
Therefore an observer moving along a line of constant (p will see all the mass of the black 
hole radiated away if — > as x~^ increases. 



3. Numerical Results For Original CGHS Model 

It may readily be seen from the form of ( |2.7| ) that for fixed x"*" they reduce to ordinary 
differential equations for the variables and dj^cf). These may then be integrated from 
large negative x~ in towards the singularity at He^'^ = 1. With and dj^cf) now known 
the solution may be evolved a step in the x"*" direction. This is of course a standard method 
in the solution of hyperbolic differential equations known as the method of characteristics. 
Here the characteristic lines are simply the x'^ and x~ directions. Details of the precise 
numerical method used are given in the appendix. 

The boundary conditions imposed are that p and </> match onto the linear dilaton 
vacuum along the line In order to specify a unique solution boundary conditions 

must also be imposed at past null infinity. In practice, the conditions are that dj^cj) and 
d^p correspond to a finite mass classical black hole solution for > Xq and at some large 
negative value of x~ where the quantum effects are always negligible. Note that A may be 
scaled out of the problem so in the following A = 1. 

To check that the numerical algorithm was functioning correctly the classical shock 
solution was evolved and the position of the apparent horizon = was plotted in 
comparison to the exact value x~ = —M. This is shown in fig. |l| where M = 50. Excellent 
agreement is obtained. 

The typical result of evolving the semiclassical equations for a large N evaporating 
solution is shown in fig. |2| . The parameters M = 50 and = 480 are thought to be 
indicative of the generic case. The apparent horizon is seen to recede as the singularity 
approaches and then crosses it. The numerical integration breaks down near the singularity 
so the path of the singularity is represented by a line of large constant R. This line begins 
to recede to larger values oi x~ as a;"*" increases signalling the appearance of a timelike 
naked singularity at x~ = x~ . The true path of the singularity lies at slightly larger values 
of x~ than this line and becomes tangent to x~ = x~ at the moment when the singularity 
becomes naked. After this point both the apparent horizon and the singularity recede 
beyond the region determined by propagation of the equations. If one follows a line of 
constant x~ close to the line x~ = x~ the curvature will increase to some maximum value 
and then decrease. As this line approaches x~ the ma:ximum value of the curvature seems 
to increase without bound. 

As expected the curvature on the apparent horizon seems to increase without bound 
as shown in fig. ^ . In fig. |^ the curvature measured by an observer who travels along a line 
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of constant (j) is shown. Equation (|2.9D relates this curvature to the effective mass of the 
black hole. This passes through zero becoming negative as the naked singularity continues 
radiating away to negative mass. 



These results are qualitatively the same as the results of [p!0H l2| in that a naked 
singularity appears. The semiclassical equations break down just before this singularity 
appears, by which time the mass of the black hole is of order the Planck mass. 

4. Another Model: Reduction of the Einstein Action From 4d to 2d 

The starting point for this model is the four dimensional Einstein action 

In J 

Restricting to spherically symmetric field configurations allows one to write the four di- 
mensional metric as 

(4)c^s2 _ ^e''Pdx+dx- + e-^'^dn^ , (4.2) 

where the dilaton field (p iiow becomes part of the four dimensional metric. The four 
dimensional Ricci curvature scalar is then 

(^^i? = 26^'^ + e-^P{8d+d-p + 2Ad+(t)d-(t) - IQd+d-cf)) . (4.3) 

Substituting this into ( [4.1| ) then gives the two dimensional dilaton gravity action 

ms=l.j d'a^/^ (^(2)^ ^ 2(V(/>)2 + 2e^^^ . (4.4) 



This model has been considered at the classical level before, for example in [|T3 . 

The static solutions of this classical action are of course the familiar Schwarzchild 
solutions. These are most easily investigated not in conformal gauge, but in the gauge 

i^)ds^ = -h(r)dt^ + T^dr^ . (4.5) 

h[r) 

Solving the equations of motion in this gauge, the most general solution up to translations 
in r obeying h{r) ^ 1 as r ^ oo is 

4>{r) =(po- log(r) 

2M (4-6) 

h{r) = 1 - — , 
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as expected. These solutions are periodic in imaginary time with period SttM, so have 
Hawking temperature T = l/8nM as in the four dimensional case. 

If one chooses to couple this classical dilaton gravity model to matter which is generally 
covariant only in two dimensions the fluctuations in the matter fields may be integrated 
out in the same way as the CGHS model to give the effective action 



Seff = ^ j dx+dx-{e-^'^[-2d+d-p + 2d+i 



2^ 



^ N 
lJ2d+fd-f + -d+pd.p 



(4.7) 



where conformal gauge has been used 

(2)c^s2 _ -e'^Pdx+dx- . (4.8) 

Ideally, one would like to include matter fields generally covariant in four dimensions. 
Integrating out the fluctuations of the matter fields in four dimensions then leads to a trace 
anomaly of the stress energy tensor that includes type terms. The effective action must 
then include highly nonlocal terms and the problem seems to become intractable. In the 
following sections the effective action ( |4.7| ) will be studied as a two dimensional toy model 
of black hole evaporation with properties much closer to four dimensional black hole physics 
than the original CGHS model. 



The equations of motion that follow from (|4.7| ) are 



24^ 



d+d-p = d+d-(f) — d+(pd-(p 
d+d.f = . 

In addition there are the constraint equations 



(4.9) 



= T++ = e-^^{4d+pd+<P - 2d\c^ + 2(a+(/>)2) + \{d^ff 

= T__ = e-2'^(4a_pa_(/)-2a2(/, + 2(a_(^)2) + i(a_/)^ 
-^((a_p)2-aip + t_(x-)) , 
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(4.10) 



where t+ and t_ are integration functions to be fixed by the boundary conditions. They 
arise from the nonlocahty stiU present in the covariant form of the effective action. 

As in the CGHS model, ([4.9| ) degenerate at a critical value of(j) = (per = — \ log(A^/24), 
where in general a curvature singularity will appear. Following similar arguments as lead 



to ( |2.8|) one finds the large approximation breaks down in this region where quantum 
gravity effects become strong. 



5. Static Solutions of the New Model 

If one is interested in finite mass static solutions it is natural to choose as the radial 
coordinate a = ^{x~^ — In these coordinates the event horizon will turn out to be at 
a = — oo. The equations of motion in the static limit become 

^"-(^T + ^^3t^((*T+e--)) ^^^^ 

p" = r - (4''? , 

and the constraint equations when fi = become 

e-2^(pV - 1<P" + urn - ^(ipr -p" + t)=0. (5.2) 

The vacuum corresponds to the solution 

<P = -\ogia), p = 0. (5.3) 

Solving the linearized equations of motion about the vacuum leads to infinite mass solutions 
when t 7^ 0, and when t = finite mass solutions with asymptotic form 

2M M 

(t) = - log(a) + log(cT) , p = . (5.4) 

a a 



The results of integrating ( |5.1| ) with boundary conditions given by (|5.4|) are shown in 
fig. H and fig. H . These are very similar to the static solutions of the CGHS model 
The coupling increases as a decreases before bouncing off the region of critical coupling 
4>cr = —\ log(A^/24), and then decreases to zero coupling as a — oo. The curvature 
R = —2e~'^'^p" rises to a maximum near this bounce region before running off to — oo 
as a ^ — oo. As M 0, (pmax <Pcn and the solution closely approaches the vacuum 
outside a region close to the origin. 
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Beyond the bounce region as a ^ — oo the exp(2(p+^)) term in ( |5.1| ) may be neglected 
allowing the equations to be solved analytically 



e ^'^ = -—p + aa + b 

J. ^ 



(5.5) 



24 24 24 
where a, 6 and c are constants. The asymptotic form of these solutions as cr ^ —oo is 



; ^'^ aa + ^ log(-ao-) 

„4:8aa/N 

V—aa 



(5.6) 



It is clear from the form of the metric that cr = —oo is an event horizon at finite proper 
distance and the curvature becomes singular there. 

Another interesting class of static solutions are those with regular event horizons. As 
expected these will correspond to infinite mass black holes supported by an incoming flux 
of radiation. In this case the most convenient radial coordinate seems to be = —x~^x~, 
which ensures that the horizon is at r = 0. The equations of motion in these coordinates 
become 



(5.7) 



and the constraint equations become 



4pV' - 2<P" + 1<P' + 2i<pr = ^e^^ ((pT -p" + Ip' + ^)- (5-8) 

The boundary conditions for a regular horizon are therefore 

<P'{0) = 0, p'(0)=0, i=0, p(0) = 0, m=<f>h, (5.9) 

where p(0) has been set to zero by a scale transformation. The solutions are parameterized 
by the value of cp on the horizon. Only when (ph < 4>cr do the solutions approach the vacuum 
( ^.3|) as r — > oo. The coupling e^*^ for a typical solution is shown in fig. |^ . 

As (ph (per the solutions approach a limiting form as can be seen from fig. |^ where 
e"'' is plotted for various values of (ph- These solutions match onto the vacuum solutions 
of (|5.7| ) as r ^ oo which are 

e""^ = a — blogr 

(5.10) 

P = log(6/r) , 
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where a and b are constants. The value of b for the hmiting form is bo ~ 6.3. 

The evaporation of a large mass black hole may be well approximated by a succession 
of these static solutions with regular horizons. As the mass is radiated away (ph will 
increase towards the critical value (per- If one expects no naked singularity to appear then 
this limiting form should approach zero temperature. In fact, this is not true, as may be 
seen by a calculation of the Hawking flux. In order that the future horizon be regular one 
must demand that t-{x~) = 0. This is analogous to the condition t = for the static 
regular horizon solution. Transforming (|5.10|) to asymptotically Minkowskian coordinates 
and a~ via 

x+ = e"^l^\ X- = -e-'"~/^o , (5.11) 



leads to 



t-{a-) = -{Ux-)+D^^.ix-)) 



N 



(5.12) 



where 

is the Schwarzian derivative and a' denotes da/dx. Thus the limiting form is unstable with 
a flnite outgoing flux. This strongly suggests that a naked singularity will form. In the 
following section this will be conflrmed by direct numerical calculation for the shockwave 
geometry. 



6. Black Hole Evaporation 

In this section the equations ( [4.9| ) will be solved numerically for incoming matter shock 
waves which form black holes, and then subsequently evaporate. The same algorithm is 
used as in Section 3. 

Now the boundary conditions are that p and (p match onto the vacuum ( |5.3D along 
the null line x~^ = 1. Again at past null inflnity (or more practically large negative x~) the 
boundary condition is that d+cp and d+p correspond to a classical black hole for x~^ > 1. 
Translated into conformal gauge, this condition is 

d+<P - + ^ ^ l-log(i(x+-x-)) , d+p^— ^. (6.1) 

x^ — X (x+ — X )^ \ J [x^ — X y 
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The boundary condition is imposed in a region where the quantum effects are always 
neghgible so one may be sure there is no extra incoming energy flux other than the finite 
energy of the initial shock wave. 

From the usual adiabatic analysis the lifetime of these black holes t ~ M^. This is 
observed in these numerical studies as a strong dependence of the lifetime on the mass, in 
contrast to the CGHS case where t ~ M. Since one is interested in the endpoint of the 
evaporation this effectively restricts one to the study of initial masses of order the Planck 
mass. 

In fig. 1^ the paths of the singularity and the apparent horizon d+(p = are displayed. 
The set of parameters chosen (M = 3 and = 2400) are believed to represent the generic 
case. These collide after a finite proper time signaling the formation of a naked singularity. 
The path of the singularity is represented by a line of large constant R. The true line of 
singularity lies at slightly larger x~ , and becomes tangent to a line of constant x~ at the 
moment when the singularity becomes naked. 



In fig. the curvature on the apparent horizon is shown, which appears to become 



singular for finite a;"*". The curvature decreases along a line of constant cp as the mass of 



the black hole is radiated away, as shown in fig. 11 . 



7. Conclusions 

In this paper, black hole evaporation in the context of two quite different semiclassical 
two dimensional dilaton gravity models has been studied. The second of these is obtained 
by the dimensional reduction of the four dimensional Einstein equations so is expected to 
closely reflect physics in four dimensions. The formation of a naked singularity appears to 
be a rather generic feature of these semiclassical calculations. This is in agreement with 
the results of the exactly solved modifled models pOHT^. However, this seems to disagree 



with recent results of Hawking and Stewart [|I4[ who, in the context of the CGHS-type 
models, claim the apparent horizon persists, and a thunderbolt singularity forms. This 
kind of endpoint is similar to a stable remnant, but the singularity spreads out to inflnity 
rather than remaining in a bounded region. This disagreement may arise because in [0 
the boundary condition on a line of constant a;~ is imposed in a region where quantum 
effects are likely to be strong. One might expect this would be equivalent to sending in 
some non-zero energy flux over a long period of time balancing the outgoing flux and 
preventing the naked singularity from appearing. In this paper the boundary condition is 
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imposed in a region where quantum effects are negligible so one can be certain the matter 
shock wave carries only a finite total energy. 

In any case, the semiclassical approximation only holds until the remaining mass of 
the black hole is of order the Planck mass. One might hope that higher order quantum 
corrections will tame the naked singularity and the residual energy might then be released 
as some kind of gamma ray burst (a rather scaled down version of the big bang naked 
singularity) . Due to the energy scale involved it is unlikely this could be detected unless it 
occurred in our close proximity. If one takes the emergence of a naked singularity seriously 
it is natural to ask whether naked singularities appear as nonperturbative solutions of string 
theory. If this is true it may make string theory unpredictive and require the introduction 
of a density matrix description as in P]. 
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Appendix A. Numerical Algorithm 

The nonhnear partial differential equations to be solved take the form 

d+d-u — f{u, d-u, d+u) (A.l) 

with boundary conditions imposed along the lines and X = Xq . For fixed x"*" 

( |A.1|) may be regarded as a first order ordinary differential equation for d^u. A fourth 
order Runge-Kutta routine with adaptive stepsize in the x~ direction is used to compute 
dj^u. The efficiency of this method relies on u being known for arbitrary x~ . This is 
accomplished using a fourth order rational function interpolation routine. The additional 
error introduced by this interpolation is more than compensated by the improved efficiency 
of the adaptive stepsize method. The derivatives d-U are computed using a second order 
backward difference method. 

Having obtained d+u as a function of x~ again a fourth order Runge-Kutta method 
with adaptive stepsize is used to evolve the solution a step in the x~^ direction. Repeating 
this procedure allows to the solution to be evolved out to large x~^. 
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Figure Captions 

Fig. 1. Position of the apparent horizon for a classical shockwave solution in the original 
CGHS model. The dashed line is the numerical result, the solid line is the exact 
result. 

Fig. 2. Position of the apparent horizon and the singularity for the original CGHS model. 

The initial mass of the black hole is M = 50 and N = 480. The apparent horizon 
recedes and eventually crosses the singularity. 

Fig. 3. Curvature on the apparent horizon for the original CGHS model, M = 50 and 
N = 480. 

Fig. 4. Curvature along a line of constant 4> ^or the original CGHS model, M = 50 and 
N = 480. 

Fig. 5. Coupling ■^e'^'^ versus the radial coordinate for static solutions of the model 
obtained from the 4d Einstein equations with mass M = 8, 10, 15, and N = 480. 

Fig. 6. Curvature versus the radial coordinate for static solutions of the model obtained 
from the 4d Einstein equations with mass M = 8, 10, 15, and N = 480. 

Fig. 7. The coupling ^e'^'^ versus the radial coordinate for a regular horizon static solu- 
tion of the model obtained from the 4d Einstein equations. 

Fig. 8. versus the radial coordinate for a regular horizon static solution of the model 
obtained from the 4d Einstein equations with 0^ = —2.67, —1.70, —1.52. Here 

(j)cr = —1.50. 

Fig. 9. Paths of the singularity and the apparent horizon for the model obtained from 
the 4d Einstein equations, with M = 3 and N = 2400. The apparent horizon 
recedes and eventually crosses the singularity. 

Fig. 10. Curvature on the apparent horizon for the model obtained from the 4d Einstein 
equations, with M = 3 and N = 2400. 

Fig. 11. Curvature along a line of constant for the model obtained from the 4d Einstein 
equations, with M = 3 and N = 2400. 
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